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A class of extrasolar giant planets - the so-called ‘hot Jupiters’^ - orhit within 0.05 AU of 
their primary stars. These planets should he hot and so emit detectable infrared radiation^. 
The planet HD 209458h^’^ is an ideal candidate for the detection and characterization of this 
infrared light because it is eclipsed by the star. This planet has an anomalously large radius 
(1.35 times that of Jupiter^), which may be the result of ongoing tidal dissipation®, but this 
explanation requires a non-zero orbital eccentricity (~0.03)®’’, maintained by interaction 
with a hypothetical second planet. Here we report detection of infrared (24 pm) radiation 
from HD 209458b, by observing the decrement in flux during secondary eclipse, when the 
planet passes behind the star. The planet’s 24 pm flux is 55±10 pjy (la), with a brightness 
temperature of 1130±150 Kelvins, confirming the predicted heating by stellar irradiation^ *. 
The secondary eclipse occurs at the midpoint between transits of the planet in front of the 
star (to within ±7 min, la), which means that a dynamically significant orbital eccentricity 
is unlikely. 

Operating cryogenically in a thermally stable spaee environment, the Spitzer Spaee Teleseope^ 
has suffieient sensitivity to detect hot Jupiters at their predicted infrared flux levels*. We observed 
the secondary eclipse (hereafter referred to as ‘the eclipse’) of HD 209458b with the 24 pm channel 
of the Multiband Imaging Photometer for Spitzer (MIPS)^°. Our photometric time series obser¬ 
vations began on 6 December 2004 at 21:29 UTC, and ended at approximately 03:23 UTC on 
7 December 2004 (5 h 54 min duration). We analyze 1696 of the 1728 10-s exposures so acquired, 
rejecting 32 images having obvious flaws. The Supplementary Information contains a sample im¬ 
age, together with information on the noise properties of the data. 
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We first verify that circumstellar dust does not eontribute signifieantly to the stellar flux. 
Summing eaeh stellar image over a 13 x 13 pixel synthetie aperture (33 x 33 aresee), we multiply 
the average sum by 1.15 to aeeount for the far wings of the point spread funetion (PSF)^\ deriving 
a flux of 21.17±0.11 mJy. The temperature of the star is elose to 6000K^^. At a distanee of 47 
pe^^, a model atmosphere''^ prediets a flux of 22 mJy, agreeing with our observed flux to within an 
estimated ~ 2 mJy error in absolute ealibration. We eonelude that the observed flux is dominated 
by photospherie emission, in agreement with a large Spitzer study of planet-bearing stars at this 
wavelength". 

Our time series analysis is optimized for high relative preeision. We extraet the intensity of 
the star from eaeh image using optimal photometry with a spatial weighting funetion'^. Seleeting 
the Tiny Tim'^ synthetie MIPS PSF for a 5000 Kelvin source at 24 /rm, we spline-interpolate it 
to 0.01 pixel spacing, rebin it to the data resolution, and center it on the stellar image. The best 
centering is judged by a least-squares fit to the star, fitting to within the noise level. The best- 
centered PSF becomes the weighting function in deriving the stellar photometric intensity. We 
subtract the average background over each image before applying the weights. MIPS data includes 
per-pixel error estimates'^, which we use in the optimal photometry and to compute errors for 
each photometric point. The optimal algorithm'^ predicts the signal-to-noise ratio (SNR) for each 
photometric point, and these average to 119. Our data are divided into 14 blocks, defined by pre¬ 
determined raster positions of the star on the detector. To check our SNR, we compute the internal 
scatter within each block. This gives SNR in the range from 95 to 120 (averaging 111), in excellent 
agreement with the optimal algorithm. For each point we use the most conservative possible error: 
either the scatter within that block or the algorithm estimate, whichever is greater. We search for 
correlations between the photometric intensities and small fluctuations in stellar position, but find 
none. We also perform simple aperture photometry on the images, and this independent procedure 
confirms our results, but with 60% greater errors. 

The performance of MIPS at 24 /rm is known to be excellent'^. Only one instrument quirk 
affects our photometry. The MIPS observing sequence obtains periodic bias images, which reset 
the detector. Images following resets have lower overall intensities (by ~ 0.1-1%), which recover 
in later images. The change is common to all pixels in the detector, and we remove it by dividing 
the stellar intensities by the average zodiacal background in each image. We thereby remove 
variations in instrument/detector response, both known and unknown. The best available zodiacal 
model'^ predicts a background increase of 0.18% during the ~ 6 h of our photometry. Because the 
star will not share this increase, we remove a 0.18% linear baseline from the stellar photometry. 
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Note that the eclipse involves both a decrease and increase in flux, and its detection is insensitive 
to monotonic linear baseline effects. 

Reliably detecting weak signals requires investigating the nature of the errors. We find that 
shot noise in the zodiacal background is the dominant source of error; systematic effects are unde¬ 
tectable after normalizing any individual pixel to the total zodiacal background. All of our results 
are based on analysis of the 1696 individual photometric measurements versus heliocentric phase 
from a recent ephemeris^® (Fig. la). We propagate the individual errors (not shown on Fig. la) 
through a transit curve fit to calculate the error on the eclipse depth. Because about half of the 
1696 points are out of eclipse, and half are in eclipse, and the SNR ~ 111 per point, the error on 
the eclipse depth should be ~ (0.009) (2°-^)/(848°'^) = 0.044% of the stellar continuum. Model 
atmospheres for hot Jupiters^’^’^^^^"^ predict eclipse depths in the range from 0.2-0.4% of the stellar 
continuum, so we anticipate a detection of 4-9cr significance. The eclipse is difficult to discern by 
eye on Fig. la, because the observed depth (0.26%) is a factor of 4 below the scatter of individual 
points. We use the known period (3.524 days) and radii^ to fit an eclipse curve to the Fig. la data, 
varying only the eclipse depth, and constraining the central phase to 0.5. This fit detects the eclipse 
at a depth of 0.26% ± 0.046%, with a reduced of 0.963, denoting a good fit. Note that the 5.6cr 
significance applies to the aggregate result, not to individual points. The eclipse is more readily 
seen by eye on Fig. lb, which presents binned data and the best-fit eclipse curve. The data are 
divided into many bins, so the aggregate 5.6a significance is much less for a single bin (SNR ~ 1 
per point). Nevertheless, the dip in flux due to the eclipse is apparent, and the observed duration is 
approximately as expected. As a check, we use a control photometric sequence (Fig. lb) to elimi¬ 
nate false positive detection of the eclipse due to instrument effects. We also plot the distribution of 
points in intensity for both the in-eclipse and out-of-eclipse phase intervals (Fig. Ic). This shows 
that the entire distribution shifts as expected with the eclipse, providing additional discrimination 
against a false positive detection. 

We further illustrate the reality of the eclipse on Fig. 2. Now shifting the eclipse curve 
in phase, we find the best-fitting amplitude and at each shift. This determines the best-fit 
central phase for the eclipse, and also further illustrates the statistical significance of the result. 
The thick line in Fig. 2a shows that the maximum amplitude (0.26%) is obtained at exactly phase 
0.5 (which is also the minimum of x^)- Further, we plot the eclipse ‘amplitude’ versus central 
phase using 100 sets of synthetic data, consisting of gaussian noise with dispersion matching the 
real data, but without an eclipse. The amplitude (0.26%) of the eclipse in the real data stands 
well above the statistical fluctuations in the synthetic data. Figure 2b shows confidence intervals 
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on the amplitude and central phase, based on the values. The phase shift of the eclipse is 
quite sensitive to eccentricity (e) and is given^^ as 5t = 2Pe cos(a;)/ 7 r, where P is the orbital 
period, and uj is the longitude of periastron. The Doppler data alone give e = 0.027 ± 0.015 
(Laughlin, G., personal communication), and allow a phase shift as large as ±0.017 (87 min). 
We find the eclipse centred at phase 0.5, and we checked the precision using a bootstrap Monte 
Carlo procedure^®. The la phase error from this method is 0.0015 (~ 7 min), consistent with 
Fig. 2b. A dynamically significant eccentricity, e ~ 0.03^’’, constrained by our 3a limit of 
6t < 21 min, requires | (u; — 7 r/ 2 ) | < 12 degrees, and is therefore only possible in the unlikely case 
that our viewing angle is closely parallel to the major axis of the orbit. A circular orbit rules out 
a promising explanation for the planet’s anomalously large radius: tidal dissipation as an interior 
energy source to slow down planetary evolution and contraction^. Because the dynamical time for 
tidal decay to a circular orbit is short, this scenario posited the presence of a perturbing second 
planet in the system to continually force the eccentricity - a planet that is no longer necessary 
with a circular orbit for HD 209458b. The infrared flux from the planet follows directly from our 
measured stellar flux (21.2 mJy) and the eclipse depth (0.26%), giving 55 ± 10 /rJy. The error is 
dominated by uncertainty in the eclipse depth. Using the planet’s known radius^ and distance*^, 
we obtain a brightness temperature T 24 =l 130 ± 150K, confirming heating by stellar irradiation^. 
Nevertheless, T 24 could differ significantly from the temperature of the equivalent blackbody (Teg), 
that is, one whose bolometric flux is the same as the planet. Without measurements at shorter 
wavelengths, a model atmosphere must be used to estimate Teg from the 24 /im flux. One such 
model is shown in Fig. 3, having Teg =1700K. This temperature is much higher than T 24 (1130K) 
due to strong, continuous H 2 O vapor absorption at 24 /rm. The bulk of the planetary thermal 
emission derives ultimately from re-radiated stellar irradiation, and is emitted at 1-4 /rm, between 
H 2 O bands. However, our 24 /im flux error admits a range of models, including some with a 
significantly lower Teg (for example, but not limited to, models with reflective clouds or less H 2 O 
vapor). 

Shortly after submission of this Letter, we became aware of a similar detection for the TrES-1 
transiting planet system^^ using Spitzer’s Infrared Array Camera^^. Together, these Spitzer results 
represent the first measurement of radiation from extrasolar planets. Additional Spitzer observa¬ 
tions should rapidly narrow the range of acceptable models, and reveal the atmospheric structure, 
composition, and other characteristics of close-in extrasolar giant planets. 
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Figure 1: Observations showing our deteetion of the seeondary eelipse in HD 209458. a, Relative 
intensities versus helioeentrie phase (seale at top) for all 1696 data points. The phase is eorreeted 
for light travel time at the orbital position of the teleseope. Error bars are suppressed for elarity. The 
gap in the data near phase 0.497 is due to a pause for teleseope overhead aetivity. The seeondary 
eelipse is present, but is a faetor of ~ 4 below the ~ 1% noise level of a single measurement, b, 
Intensities from a, averaged in bins of phase width 0.001 (seale at top), with la error bars eomputed 
by statistieal eombination from the errors of individual points. The red line is the best-fit seeondary 
eelipse eurve (depth=0.26%), eonstrained to a eentral phase of 0.5. The points in blue are a eontrol 
sequenee, summing intensities over a 10 x 10 pixel region of the deteetor, to beat down the random 
errors and reveal any possible systematie effeets. The eontrol sequenee uses the same deteetor 
pixels, on average, as those where the star resides, but is sampled out of phase with the variations 
in the star’s raster motion during the MIPS photometry eyele. c. Histograms of intensity (lower 
abseissa seale) for the points in a, with bin width 0.1%, shown separately for the out-of-eelipse 
(blaek) and in-eelipse intervals (red). 
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Figure 2: Amplitude of the secondary eclipse versus assumed central phase, with confidence 
intervals for both, a, The darkest line shows the amplitude of the best-fit eclipse curve versus the 
assumed central phase (scale at top). The overplotted point marks the fit having smallest which 
also has the greatest eclipse amplitude. The numerous thinner black lines show the effect of fitting 
to 100 synthetic data sets containing no eclipse, but having the same per-point errors as the real 
data. Their fluctuations in retrieved amplitude versus phase are indicative of the error in eclipse 
amplitude, and are consistent with a =0.046%. Note that the eclipse amplitude found in the real 
data (0.26%) stands well above the error envelope at phase 0.5. b. Confidence intervals at the 
1, 2, 3 and 4a levels for the eclipse amplitude and central phase (note expansion of phase scale, 
at bottom). The plotted point marks the best fit (minimum x^) with eclipse depth of 0.26%, and 
central phase indistinguishable from 0.5. The center of the eclipse occurs in our data at Julian day 
2453346.5278. 


Figure 3: Flux from a model atmosphere shown in comparison to our measured infrared flux at 
24 /im. A theoretical spectrum (solid line) shows that planetary emission (dominated by absorbed 
and re-radiated stellar radiation) should be very different from a blackbody. Hence, models are 
required to interpret the 24 /im flux measurement in terms of the planetary temperature. The 
model shown has an equivalent blackbody temperature =1700 K and was computed from a 
one-dimensional plane-parallel radiative transfer model, considering a solar system abundance of 
gases, no clouds, and the absorbed stellar radiation re-emitted on the day side only. Note the 
marked difference from a 1700K blackbody (dashed line), although the total flux integrated over 
the blackbody spectrum is equal to the total flux integrated over the model spectrum. (The peaks 
at short wavelength dominate the flux integral in the atmosphere model, note log scale in the 
ordinate.) The suppressed flux at 24 /im is due to water vapor opacity. This model lies at the hot 
end of the range of plausible models consistent with our measurement, but the error bars admit 
models with cooler T^q. 
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